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A new analysis is given of numerical simulation data on the archetype square lattice Ising Spin 
Glasses (ISG) with a bimodal (± J) and Gaussian interaction distributions. It is well established that 
the ordering temperature of both models is zero. The Gaussian has a non-degenerate ground state so 
exponent 77 = 0 and it has a continuous distribution of energy levels. For the bimodal model, above 
a size-dependent cross-over temperature T*{L) there is a regime of effectively continuous energy 
levels; below T* (L) there is a distinct regime dominated by the highly degenerate ground state plus 
an energy gap to the excited states. T*{L) tends to zero at very large L leaving only the effectively 
continuous regime in the thermodynamic limit. We show that in this regime the critical exponent 
rj is not zero, so the effectively continuous regime 2D bimodal ISG is not in the same universality 
class as the 2D Gaussian ISG. The simulation data on both models are analyzed using a scaling 
variable r = r^/(l -|- T^) suitable for zero temperature transition ISGs, together with appropriate 
scaling expressions. Accurate simulation estimates can be obtained for the temperature dependence 
of the thermodynamic limit reduced susceptibility x(''') second moment correlation length ^(r) 
over the entire range of temperature from zero to infinity. The Gaussian critical exponent from the 
simulations u = 3.5(1) is in full agreement with the well established value from the literature. The 
bimodal exponent from the thermodynamic limit regime analysis is = 4.2(1), once again different 
from the Gaussian value. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 


I. INTRODUCTION 

The canonical Edwards-Anderson (EA) model Ising 
Spin Glasses (ISGs) in dimension d = 2 have been the 
subject of very many numerical studies. There is now 
consensus supported by analytic arguments that the two 
archetype models, the ISGs on square lattices with near 
neighbor interactions having distributions which are ei¬ 
ther Gaussian or Bimodal (±J), have zero temperature 
transitions Qi- For the Gaussian case, where the in¬ 
teraction distribution is continuous and the ground state 
is unique, there is now also general consensus concern¬ 
ing the low temperature thermodynamic limit (ThL) be¬ 
havior and exponents. In the bimodal case there is an 
’’effectively continuous energy level distribution” regime 
coming down from high temperatures and ending with 
a crossover at a size dependent temperature T*{L) to a 
ground state dominated regime Q. Interpretations dif¬ 
fer considerably concerning the critical exponents for the 
bimodal interaction case. We will give a new analysis of 
accurate numerical Monte Garlo data up to size L = 128 
on the bimodal system in the ThL and up to L = 64 on 
the Gaussian system. We use the temperature T or the 
inverse temperature /3 when it is convenient. 

We first discuss discuss the specific heat using data 
from the simulations and independent data down to low 
temperatures and large sizes from Refs. Q and Q. Then 
we analyse the simulation data for other observables to 
obtain reliable and accurate estimates for the critical ex¬ 
ponents of the ThL regime, using both the conventional 
scaling variable T and a novel scaling variable compatible 
with the generic scaling approach for ISGs introduced in 
0, adapted to a situation where = 0. We find that 


the Gaussian model and the bimodal model in the ThL 
are not in the same universality class. 

The 2D Gaussian model is relatively clear-cut. Because 
Tc = 0 and the interaction distribution is continuous, 
there is a unique ground state (for each sample) and the 
low temperature excitation distribution has no gap. The 
fact that the ground state is unique necessarily implies 
that for all A, as r —>• 0, (T, A) —)• 00 and the reduced 
susceptibility x(T, A) —)• L?. With T chosen as the crit¬ 
ical scaling variable, the standard thermodynamic limit 
low temperature critical expressions are ^(T) ~ and 
x{T) ~ T~'^ = because the critical exponent 77 

is strictly zero. The critical behavior of both the ob¬ 
servables at low temperature is governed by the single 
exponent v, which is related to the stiffness exponent 
through 6 = —Xjv. Accurate zero temperature domain 
wall stiffness measurements to large sizes [3, [Mil show 
that B = -0.285(2), i.e. = 3.50(3). 

In the 2D bimodal case the situation is complicated by 
two factors. First, the ground state is not unique but is 
massively degenerate; the zero-temperature entropy per 
spin is So = 0.078(5)^5 [llHl^l- Secondly, the distribu¬ 
tion of excited state energy levels is not continuous but 
increases by steps of 4J; in particular there is an energy 
gap 4J between the ground state and the first excited 
state. One can write Q the ’’naive” leading low temper¬ 
ature finite size specific heat expression 

n fa T\ _ 16J2exp(S'i(A)-S'o(i))exp(-4J/r) 

Cv[p,L) — 

where S'i(A), S'o(A) are the sample-averaged entropies 
of the first excited state and the ground state. Setting 
J = I a crossover temperature can then be defined by 
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I,Hill 

T*{L) = A/iS,{L)-So{L)) (2) 

which separates the critical behavior in the low temper¬ 
ature ground state dominated regime (with Cv{/3,L) ^ 
exp(—4/T)/T^) andaT > T*{L) regime where the whole 
ensemble of higher energy states dominate the thermody¬ 
namics Q. An explicit phenomenological expression for 
T*{L) derived from Eq. 5 of Ref. [l^, which is consistent 
with the raw data points [1,[I^ for Si{L) — So{L) is 


exp(0.199 ^ -I- L‘^{\n{L‘^) — 1)) -|- 0.473) 

A much simpler droplet-based expression from Q is 
T*{L) ss T*[L) decreases with increasing L be¬ 

cause the degeneracy of the excited states increases faster 
with L than that of the ground state. We will assume Q 
that in the T > T* (L) ThL regime the data can be anal¬ 
ysed in the same way as if the energy level distribution 
were continuous. With this assumption the T > T*{L) 
regime will have ” effectively continuous” energy level dis¬ 
tribution critical exponents with an effective ordering 
temperature Tc still zero. The ground state dominated 
regime at T < T*{L) is a finite size effect which disap¬ 
pears in the infinite L limit. 

A non-zero r] is to be expected a priori for a system 
with a strong ground state degeneracy, unless each indi¬ 
vidual ground state is isolated in phase space which is 
not the case m in the bimodal ISG. A droplet anal¬ 
ysis of g round state measurements on large sized sam¬ 
ples [Ig show that rj ss 0.22, broadly consistent with a 
number of finite temperature simulation estimates [Ti¬ 
ll^. However it has been claimed that in the T > T*{L) 
regime the bimodal ISG can be considered to be effec¬ 
tively in the same universality class as the Gaussian ISG 
Q , meaning that the effective exponents are again r/ = 0 
and r' = 3.50(3). In view of the basic definition of rj in 
terms of the short range limit of the spin-spin correla¬ 
tion function G(r, T) = G[r~^ exp(r/f (T))], this claim is 
rather surprising. 

A major difficulty in establishing the limiting [T > 
T*(L),L ^ oo,T —0] behavior for ISGs in dimension 
2 |14| consists in finding an appropriate and reliable ex¬ 
trapolation procedure from simulation data necessarily 
restricted in size and in temperature because of the need 
to achieve good thermal equilibration at large sizes. This 
is a problem that we will address. 

II. SIMULATIONS 

The simulations were performed using the Houdayer 
cluster method 18| in combination with the exchange 
Monte Carlo [M| method. In the cluster step we first 
pick a random site i and compute its overlap qt = SfSf, 
where Sf and denote the spin for two different repli¬ 
cas. We then build an equal-g cluster along the nearest 


neighbor interactions and flip all cluster spins in both 
replicas. We used four replicas which turned out to be 
remarkably efficacious. On each iteration the replicas are 
paired at random, then, for each pair, a cluster update 
is performed, and the usual heat-bath spin update and 
exchange. 

For all systems we used /3max = 3.0. The number 
of temperatures were more than 250 for the smallest 
systems starting at /3min = 0.2. With increasing sys¬ 
tem size the number of temperatures was decreased and 
/3min increased. For the largest system (bimodal Jy with 
L = 128) 70 temperatures were used with /3jnin = 1-2. 
The exchange rate was always at least 0.3 for all sys¬ 
tems and temperatures. The systems were deemed equi¬ 
librated when the average (g^) for the systems at /3max 
appeared stable between runs. The number of equilibra¬ 
tion steps increased with system size, for the bimodal 
L = 128 this took about 600000 steps. After equili¬ 
bration, at least 200000 measurements were made for 
each sample for all sizes, taking place after every cluster- 
sweep-exchange step. 

The usual observables were registered, the energy 
E{/3,L), the correlation length the spin over¬ 

lap moments {\q\), (q^), (|?P), {q'^)- Correlations 
{E{I3,L),U{I3,L)) between the energy and some observ¬ 
ables U{j3,L) were also registered. Thermodynamic 
derivatives could then be evaluated through the usual 
aC/(/3,L)/5/3 = (C/(/3,L), A(/3,L)) - (C/(/3, L))(A(/3, L)). 
Error estimates of observables and derivatives were done 
with the bootstrap method. 

Sizes studied were L = 4, 6, 8, 12, 16, 24, 32, 48, 
64, 96, 128 for bimodal interactions, and up to L = 64 
for the Gaussian interactions, with 2^^ = 8192 samples 
(Jij-interactions) for all sizes. 

III. SPECIFIC HEAT 

The size dependence of the ground state energy per 
spin e(0, L) for the 2D Gaussian ISG has been shown 
[l^ to follow the simple critical finite size scaling 
rule 

e( 0 ,L)-e( 0 ,oo) - ( 4 ) 

with a 9 consistent with the estimate from ground state 
domain wall stiffness measurements [I| . Standard scaling 
arguments [J would suggest that the low temperature 
specihc heat should behave as 

a(/3,L)~/3-2‘'«/3-" (5) 

but because of the continuous interaction distribution, 
in addition to critical excitations there are always sin¬ 
gle spin excitations. These lead to a term C'„(/3,L) « T 
which dominates the Gaussian low temperature specific 
heat as noted by Ref. 0. Specific heat data for the bi¬ 
modal model were calculated through the present simula¬ 
tions; data extending to a much lower temperature range 
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and larger sizes have already been measured using the 
sophisticated Pfaffian arithmetic technique by Lukic et 
al and by Thomas et a/ and we are very grateful 
to be able to quote these results in extenso. 

The data for the two models are shown, see Fig. [T] 
and Fig. [21 in the form of plots of the derivative y = 
91n(C'„(/3, L))/9/3 against x = T. This non-conventional 
form of plot happens to be particularly instructive. A low 
temperature limit Cv{P,L) ~ appears as a straight 
line through the origin with slope —x, while a low tem¬ 
perature limit of the ” naive” bimodal ground state dom¬ 
inated form Eq. © appears as a straight line with inter¬ 
cept —4 and slope +2. 

The Gaussian data are almost independent of L for 
the whole temperature range. Physically this occurs be¬ 
cause the specihc heat in ISGs is predominately a near 
neighbor effect. The curve tends to a slope dy/dx —1 
corresponding to Cy ^ in the low T limit, in agree¬ 
ment with the conclusion of Ref. Q . 

For the bimodal model there is first a high tempera¬ 
ture and/or high L envelope curve corresponding to the 
effectively continuous T > T*{L) regime. In this regime 
finite size effects are very weak : the specific heat is al¬ 
most independent of L as in the Gaussian. The curves 
for the two models are of similar form but are not iden¬ 
tical. In the large L, low T limit of this envelope curve, 
the bimodal data as shown in Fig. |2| indicate Cy ^ in 
agreement with the conclusions drawn in Ref. Q based 
on droplet excitation arguments. 

For each L the data curve peels off the large L enve¬ 
lope curve below an L-dependent temperature which can 
be identified with the start of the effectively continuous 
to ground state dominated regime crossover centered at 
T*{L). Finally for each L in the low temperature range 
T <g; T*{L) the specific heat links up to the ’’naive” 
limit of Eq. (|T|). (It should be noted that because of 
the logarithmic derivative, temperature independent in¬ 
dependent factors do not show up in this plot). The 
crossover can be seen to be gentle for small L, becom¬ 
ing sharp for large L. Defining T*{L) as the location of 
the maximum positive slope on this plot, the crossover 
temperatures can be clearly identified and are consistent 
with = 1 . 1 ( 1 ). 

An anomalous limit of the form 

C'„(/3) ~/?^exp(-2/3) ( 6 ) 

which has been proposed by some authors [1, fol¬ 
lowing Ref. [ 2 ^ is inconsistent with the data in Fig. |2| 
for all L and T (see also [13, H^). An intermediate L 
regime where Cy{f3,L) ~ 7"5-25 g^g proposed in Ref. [11] 
or Cv{l3,L) ^ as proposed in Ref. 0 appear to be 
valid only for a limited range of T and L. 

IV. THE EXPONENT y 

For ISGs with non-zero critical temperatures finite size 
scaling analyses at and close to the critical temperature 


are used to estimate critical exponents. For the 2D bi¬ 
modal ISG, because of the crossover to the ground state 
dominated regime this approach is ruled out and the crit¬ 
ical exponents must be estimated using the entirely dif¬ 
ferent strategy of ThL measurements. 

The standard Renormalization Group theory (RGT) 
scaling variable for models with non-zero ordering tem¬ 
peratures is t = (T — Tc)/Tc. This obviously cannot be 
used when Ty = 0 ; by convention the scaling variable 
used in the literature for 2D ISGs is the un-normalized 
temperature T. This is only a convention; it is perfectly 
legitimate to use other conventions. Thus, when con¬ 
sidering the canonical ID Ising ferromagnet, Baxter [llj 
remarks ’’When Tc = 0 it is more sensible to replace 
t = (T — Tc)/Tc by T = exp(—2/3)”. (In fact for the 
particular ID model scaling without corrections over the 
entire temperature range follows if a related scaling vari¬ 
able T = 1 — tanh(/3) is chosen [13, H^). Below we will 
introduce another scaling variable appropriate for ISGs 
with Tc = 0 , but for the moment we follow this tradi¬ 
tional t = T 2D ISG convention. The critical exponents 
are defined through the leading ThL expressions for the 
reduced susceptibility and the second moment correla¬ 
tion length within this convention : x{T) = 
and ^(T) = in the limit T —>• 0,L —>• 00 . For 

all data which fulfil the condition (either in the bimodal 
and Gaussian models) L > K^{T,L) with K ob¬ 

servables such as x(T, L) and ^{T,L) depend on T but 
not on L, and so correspond to the ThL infinite size val¬ 
ues x(r) and ^{T). The ThL condition defines implic¬ 
itly a crossover temperature T^{L). It turns out that in 
the bimodal 2D ISG the ThL limit temperature T^[L) 
is always higher than the corresponding crossover tem¬ 
perature to the ground state dominated regime T*{L) 
defined above, so the ThL data are always well in the 
effectively continuous regime. The ThL data extrapola¬ 
tion to T = 0 corresponds to estimates for the critical 
exponents in the successive limits [L oo,T —>■ 0] so 
in the effectively continuous energy level regime, to be 
distinguished from the exponents defined taking the suc¬ 
cessive limits [T —> 0,L —> 00 ] which would correspond 
to the ’’finite size” ground state dominated regime. 

There have been many previous studies having the aim 
of estimating the critical exponents and in particular y 
for the bimodal model. McMillan already in 1983 esti¬ 
mated 7 / = 0.28(4) from G{r) correlation data on one 
L = 96 sample well in the effectively continuous regime 
0 . Katzgraber and Lee [l9l|e stimated 77 = 0.138(5) 
from x{T,L) data. Jorg et al 3 show a plot of Inx(T) 
against In^(r) after an extrapolation to infinite L using 
the technique of Ref. [ 2 ^ . They state ” fits of this curve 
lead to values of 77 that are very small, between 0 and 
0.1, strongly suggestive of 7/ = 0”. However, this type of 
extrapolation to infinite L is delicate, particularly in the 
bimodal 2D case. 

In addition, the data displayed by Q on a Inx(r) — 
ln,f(r) plot extending over five decades on the y axis 
are hard to fit with precision. Katzgraber et al [l3| 
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show a plot of 91nx(T, L)/91n^(T, L) which in princi¬ 
ple is equivalent to the Ref. Q plot but which provides 
a display much more sensitive to the value of 77 ; they 
state cautiously ’’for all system sizes and temperatures 
studied T^eff is always greater than 0 . 2 , although an ex¬ 
trapolation to 77 = 0 cannot be ruled out”, so that the 
possibility of the bimodal and Gaussian ISGs being in 
the same universality class ” cannot be reliably proven”. 
In Refs. [H, it is claimed that the Gaussian and bi¬ 
modal models are in the same universality class, which is 
surprising as ’’the data are not sufficiently precise to pro¬ 
vide a precise determination of 77 , being consistent with 
a small value 77 < 0 . 2 , including 77 = 0 ”. 

All the estimates quoted so far can be considered to 
concern the effectively continuous regime. At zero or low 
temperatures, so in the ground state dominated regime, 
different sophisticated algorithms lead to the estimates 
77 = 0.14(1) [13, and to 77 = 0.22 [H]. 

In Fig. [3] and Fig. SI we show plots of 77 = 

i91nx(T, L)/91n^(T, A) against x = 1/^{T,L) for the 
Gaussian and bimodal models. These are raw data points 
having the high statistical precision of the present mea¬ 
surements. With the conventional definition of the criti¬ 
cal exponents through x(r, L) and ^(T, L) ^ T’-' 

in the ThL regime low-T limit, the limiting slope dy/dx 
at criticality as a; ^ 0 is by definition equal to 2 — 77 . For 
the Gaussian model the observed tendency of the slope 
is consistent with the limit of 77 = 0 which must be the 
case for this nondegenerate ground state model. For the 
bimodal model the observed y(x) in the ThL regime is 
not tending to 2 but to a constant limit of 1.80(2). Slight 
overshoots for each L in both systems can be ascribed to 
x{T,L) and ^{T,L) not reaching the ThL condition at 
quite the same temperature. As stated above, very sim¬ 
ilar observations were made in Ref. 0 for the bimodal 
model. The present results thus confirm unambiguously 
that for the bimodal ISG in the effectively continuous 
regime 77 is not zero but is w 0.20. Thus the bimodal 
ISG in the effectively continuous ThL regime and the 
Gaussian ISG are not in the same universality class. 

It has been shown that in dimension 4 also, Gaussian 
and bimodal ISGs are not in the same universality class 
[sJl, so the breakdown of universality in ISGs appears to 
be general. 


V. SCALING AND ZERO TEMPERATURE 
ORDERING 

Estimating the exponents tz or 7 = (2 — 77 ) 7 ^ is more 
difficult than for the exponent 77 . As we have noted 
above, the standard RGT convention for models with 
finite temperature ordering is to use the scaling vari¬ 
able t = (T/Tc) — 1, which obviously cannot be ap¬ 
plied to models with Tc = 0, and for 2D ISGs the pre¬ 
ferred convention in the literature has been to use the 
un-normalized scaling variable t = T. In practice this is 
inefficient as the extrapolations towards the T = 0 limit 


in order to estimate the values for the critical exponents 
are very ambiguous. For instance, when presenting T- 
scaled susceptibility data for sizes up to A = 128 Katz- 
graber et al [l 3 | state ’’the [susceptibility] data for the 
bimodal case can be extrapolated to any arbitrary value 
including 1 / Xea = 0 ”. 

We will introduce a novel scaling variable suitable for 
the 2D ISGs, applying the same principles as for ISGs at 
higher dimensions ( 3 , adapted to A = 0 ordering : 

- For spin glasses the relevant interaction strength pa¬ 

rameter is not J but is (J^), so the natural dimensionless 
parameter is (or alternatively tanh^(J/3) for bi¬ 

modal ISGs). With the standard normalisation (J^) = 1 
the natural inverse ’’temperature” in ISGs is /3^, not 
/3. This was recognized immediately after the Edwards- 
Anderson model was introduced, in high temperature se¬ 
ries expansion (HTSE) analyses for ISGs including 2D 
models but has since been overlooked in most 

simulation analyses. 

- It is convenient to choose a scaling variable r defined 
in such a way that r = 0 at criticality and r = 1 at 
infinite temperature. With an ISG ordering at a finite 
inverse temperature Pc, t{P) = 1 — P'^/Pi is an appropri¬ 
ate choice @,[ 33 . When Pc = 00 as in the 2D ISG case, 
Tt{P) = I —(tanh(/3)/tanh(/3c))^ = I—tanh(/3)^ has been 
used [13, but here we will prefer Tb{P) = 1/(1 -1-/3^) as it 
turns out to be efficient and the limits are easy to relate 
to those of the T scaling convention. With non-zero Tc 
the effective exponents at criticality do not depend on 
the choice of scaling variable; this is not the case when 
Tc = 0, but a simple dictionary is given below relating 
the limiting derivatives for Tf, scaling to the exponents for 
the conventional T scaling. 

- The ThL HTSE Darboux (sl] format for observables 
Q{x) is 

Q{x) = 1-1- aix + 021 ^ -I- asx^ H- (7) 

with X = P^ in ISGs [33. The HTSE ISG susceptibility 
xiP'^) is naturally in this format, so for ISG models with 
Tc > 0 the ThL susceptibility can be scaled in the Wegner 
[Inform 

xiP^) = C'^t(/32)-G-'^)-F[1 + a^riP^f + ■ • ■ ] (8) 

Because the correlation function second moment 772 
HTSE is of the form (see Ref. [13 for the Ising ferro- 
magnet) 

H2{x) = X + aix^ + a2X^ + ■ ■ ■ ( 9 ) 

and the second moment correlation length is dehned 
through ^(a;)^ = y, 2 /{zx{x)) with z the number of near¬ 
est neighbors, for consistency the appropriate correlation 
length variable for ISG scaling is ^{x)/P rather than f (x) 
(whether Tc is zero or not). This point has been spelt out 
in Ref. 0. 

Examples of applications of the scaling rules outlined 
here to other specific models (both ferromagnets and 
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ISGs) have been given elsewhere. A general discussion of 
ferromagnets and spin glasses is given in Ref. Q , analyses 
of 3D Ising, XY and Heisenberg ferromagnets in Ref. , 
the 2D Ising ferromagnet is analysed in Ref. [42| , 3D Ising 
ferromagnets in , high dimension Ising ferromag¬ 

nets in Ref. [i^ , and the 2D Villain fully frustrated model 
in Ref. 

The scaling of the Binder cumulant 





( 10 ) 


is discussed in the Appendix. The 2D simulation data 
analysis and the extrapolations below are based on the 
derivatives dhiQ{Tb, L)/dlnri, in the ThL regime where 
these derivatives are independent of L and so equal to the 
infinite size derivatives. An advantage of the 2D models 
is that in contrast to t(/ 3^) for the models with non-zero 
Tc, for the 2D ISGs with Tc = 0 there is no uncertainty 
in the definition of Tf,(/?^) related to an uncertainty in the 
value of the ordering temperature. 

Once the Tf, —^ 0 limits for the various derivatives have 
been estimated by extrapolation of the ThL data for fi¬ 
nite L, there is a simple dictionary for translating into 
terms of the conventional T scaling critical exponents v 
and r] defined above : 


_ ^liix(T'b) - v) 

i9 In Tb 2 

dHTCin)) (ly-l) 

d In Tf, 2 

'91nx(Tb) -7?) 

d\n{T^{Tb)) v-\ 

^ In g(rb) 
d\nTb 


(lla) 

(llb) 

(llc) 
(lid) 


VI. ANALYSES WITH THE SCALING 
VARIABLE n 

The four derivatives of Eq. (fTTl) are shown in Figs.[S]to 
In contrast to the derivatives in which T is used as 
the scaling variable, each derivative can be extrapolated 
in a fairly unambiguous manner to criticality, and always 
has an exact finite value at infinite temperature p, = 1 . 

The exact infinite temperature limits from the general 
high temperature scaling expansion expressions (37| ap¬ 
plied to scaling with rj, are (when p, —1): 

_ 91nx(rb) _ ^ ... = 4 (bimodal) (12a) 

alnrb 

= 1 (Gauss.),... = - (bimodal) (12b) 
OlllTb 3 

ainTO?))) = ^ ^ = f 

The extrapolation method is outlined in Appendix II. 
With r] = 0 and assuming v = 3.48(5) m , the predicted 


Gaussian critical limits (when Tb —>■ 0) for the derivatives 
are 


_ dlnxin) 

dluTb 

dHmn)) 

d\nTb 

dlnxin) 

dHmn)) 

^lnff(Tb) 

dlriTb 


3.48(5) 

1.24(3) 

2.81(1) 

3.5(3) 


(13a) 

(13b) 

(13c) 

(13d) 


From the fitted ThL data extrapolations (see Appendix 
II) the estimated Gaussian critical limit values are 
3.40(10), 1.25(5), 2.90(10), 3.6(1) respectively. These val¬ 
ues are fully consistent with the list above, validating the 
approach and the extrapolation procedure that we have 
used. For the bimodal data, the extrapolated ThL limits 
(wheu Tf, —0) from the figures (see Appendix II) give 
estimates 


_ dlnxin) 

d\nTb 

dHmn)) 
d In Tb 
dlnxin) 
dHmn)) 

d\ng{Tb) 

d\nTb 


^4.3(1) 

(14a) 

^1.9(1) 

(14b) 

2.15(10) 

(14c) 

^ 4.7(1) 

(14d) 


These are all significantly different from the Gaussian val¬ 
ues, confirming non-universality. When translated into 
the T scaling convention, the 2D bimodal critical ex¬ 
ponents from these measurements are g = 0.20(2) and 
V = 4.8(2) (so 7 = (2 — g)i/ = 8.6(4)). Not only is 
the ThL bimodal exponent g different from the Gaussian 
value but the v value is different also. 


VH. CONCLUSION 

Simulation data are presented for the Gaussian and 
bimodal interaction distribution Ising spin glasses in di¬ 
mension two, which are known to order only at zero 
temperature. In order to facilitate extrapolations to 
zero temperature, a temperature scaling variable Tf, = 
T^/(l -I- T^) is introduced in addition to the traditional 
2D ISG scaling variable t = T. The Gaussian simulation 
data are completely consistent with the well established 
critical behavior for this model, with exponents rj = 0 
and v = 3.48(5) M- 

The bimodal specific heat simulation data supple¬ 
mented by data from Lukic et al Q and from Thomas et 
al Q show clear crossovers from an effectively continuous 
energy level thermodynamic limit regime to a finite size 
ground state dominated regime at size dependent tem¬ 
peratures T*{L) ss 1.1/L^/^ (see Ref. @). 





















6 



T 


FIG. 1. (Color on line) Gaussian 2D ISG. Logarithmic deriva¬ 
tive of the specific heat d\nCv{T, L)/djS against T. Sizes 
L = 64, 48, 32, 24, 16, 12 top to bottom in the dip. Curve : 
extrapolation. 



FIG. 3. (Color on line) Gaussian 2D ISG. Derivative 
d\nx{T,L)/d\n^{T,L) against 1/^{T,L) for L = 64, 48, 32, 
24 left to right. In this and all following figures both Gaus¬ 
sian and bimodal, the color coding is : black, pink, red, blue, 
green, brown, cyan, olive for L — 128, 96, 64, 48, 32, 24, 16, 
12 . 




FIG. 2. (Color on line) Bimodal 2D ISG. Logarithmic deriva¬ 
tive of the specific heat d\nCv{T, L) /df5 against T. Full 
points : simulation data L = 96, 48, 24, 12, 8 (black, green, 
red, pink, cyan) top to bottom. Open points : Pfaffian data 
; red polygons L = 512, blue right triangles L = 256, black 
left triangles L = 128, brown diamonds L = 64 (all data 
from Ref. i). green down triangles L = 50, red m triangles 
L — 24, pink circles L — 12, all data fom Ref. [j. Dashed 
diagonal red line y — —3x, green diagonal line y = —A + 2x. 


The extrapolated ThL simulation results tend to crit¬ 
ical limits which correspond consistently to rj = 0.20(2) 
and jz = 4.8(2), clearly different from the Gaussian val¬ 
ues. This demonstrates that the standard universality 
rules do not hold for 2D ISG models. In dimension 4 
also bimodal and Gaussian ISGs have been shown not to 
be in the same universality class either , strongly sug¬ 
gesting a lack of universality for ISGs in each dimension 
(presumably up to the upper critical dimension). 


FIG. 4. (Color on line) Bimodal 2D ISG. Derivative 
dlnx{T, L)/dln^{T, L) against 1/^{T,L) for L = 128, 96, 
64, 48, 32, 24 left to right. Same color coding as in Fig. |3l 


VIII. APPENDIX I: BINDER CUMULANT 


The ferromagnetic Binder cumulant has been exten¬ 
sively exploited in the FSS limit regime very close to 
criticality for its properties as a dimensionless observ¬ 
able. In addition its ThL properties can also be studied. 
In Ising ferromagnets, the critical exponent for the sec¬ 
ond field derivative of the susceptibility X4 (also called 
the non-linear susceptibility), is 

74 = 7 -b 2A = i/d + 2x (15) 


The non-linear susceptibility X 4 is directly related to the 
Binder cumulant, Eq. 10.2, through 




-X4 

Ldx^ 


3{m'^) — (m^) 

{mP'Y 


(16) 
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FIG. 5. (Color on line) Gaussian 2D ISG. The derivative 
d\nx{T, L) / dluTb against Tb. Sizes L = 64, 48, 32, 24, 16 
left to right. Same color coding as in Fig. [3] Dashed line: 
extrapolation. Red arrow : exact infinite temperature value. 
Blue arrow : Gaussian critical value. 



FIG. 6. (Color on line) Bimodal 2D ISG. The derivative 
dlnxiT, L)/d\nTb against Tb- Sizes L = 128, 96, 64, 48, 32, 
24, 16 left to right. Same color coding as in Fig. [3] Dashed 
line: extrapolation. Red arrow : exact infinite temperature 
value. 

As X scales with the critical exponent 7 , the normalized 
Binder cumulant L) scales with the ThL regime 

critical exponent d\n{L‘^g)/dlnr = {vd + 27 ) — 27 = vd. 
In any S' = 1/2 Ising system the infinite temperature 
(independent spin) limit for the Binder cumulant is 

g{^,N) = l/N (17) 

where N is the number of spins; N = L‘^ for a hyper- 
cubic lattice. So L‘^g{l3, L) has an infinite temperature 
limit which is strictly 1, and a large L critical limit (with 
corrections as for the other observables): 

+ (18) 



FIG. 7. (Color on line) Gaussian 2D ISG. The derivative 
dln{T^{T, L))/dlnTb against Tb- Sizes L = 64, 48, 32, 24 
left to right. Same color coding as in Fig. [3] Dashed line: 
extrapolation. Red arrow : exact infinite temperature value. 
Blue arrow : Gaussian critical value. 



FIG. 8. (Color on line) Bimodal 2D ISG. The derivative 
dln{T^{T, L))/dlnTb against Tb- Sizes L = 128, 96, 64, 48, 
32, 24, 16 left to right. Same color coding as in Fig. [S] Dashed 
line: extrapolation. Red arrow : exact infinite temperature 
value. 

Exactly the same argument can be transposed to ISGs 
(see Ref. [S^ for Xi in ISGs). In the particular case of a 
2D ISG model with r;, scaling, the critical value for the 
derivative d\n{L‘^g{Tb, L))/dlnTi, of the Binder cumulant 
ThL data extrapolated to Tf, = 0 is 2v(2 = v where v is 
once again the correlation length critical exponent in the 
T scaling convention. The Binder cumulant data plotted 
in the Eq. form are shown for the two models in 
Figs. [TSl and [Tdl The ThL envelope curves can be seen by 
inspection. The derivatives of these curves have already 
been shown in Figs. 191 and [TUI 

It has been suggested that if two models have the 
same function when y = g{f3, L) is plotted against x = 
^{I3,L)/L, it is a proof of universality. However, because 
both Lg^/'^[j3,L) and ^(/3,L) are controlled by just the 
same exponent u this is questionable. 


L'^gin, L) 






















FIG. 9. (Color on line) Gaussian 2D ISG. The derivative 
d\ng{T, L)/dhiTb against Tb, where g{T,L) is the Binder cu- 
mulant. L = 64, 48, 32, 24 left to right. Same color coding as 
inFig.|3] Line: extrapolation. Blue arrow : Gaussian critical 
value. 



FIG. 10. (Golor on line) Bimodal 2D ISG. The derivative 
dlng{T, L)/d\nTb against Tb, where g{T,L) is the Binder cu- 
mulant. Sizes L = 128, 96, 64, 48, 32, 24 left to right. Same 
color coding as in Fig. [3] Line: extrapolation. 


IX. APPENDIX II 

As the data sets do not extend to infinite size, to es¬ 
timate the critical = 0 limit values from the ThL 
derivative data in Figs. [5] to [121 an extrapolation must 
be made. There is no definitive method to extrapolate 
so as to be sure to obtain exact values of the critical 
exponents, though data to still larger sizes would make 
the task easier. The most economical choice for extrap¬ 
olation is to assume that the ThL derivative data con¬ 
tinue to evolve smoothly and regularly when an extrap¬ 
olation is made towards = 0 through the smaller tj, 
region where no ThL data are for the moment avail¬ 
able. To do this, for each derivative observable y{x) 



FIG. 11. (Golor on line) Gaussian 2D ISG. The derivative 
d\nx{T,L)/d\n{T^{T,L)) against Tb- Sizes L = 64, 48, 32, 
24 left to right. Same color coding as in Fig. (3) Dashed line: 
extrapolation. Red arrow : exact infinite temperature value. 
Blue arrow : Gaussian critical value. 



t. 2.5- 
€ 

2.3- 

2.2 H-1-1-1-1-1-1-1-^-1-1-1-1-1-^-1-1-1-1-1-1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


FIG. 12. (Golor on line) Bimodal 2D ISG. The derivative 
dlnx{T, L)/dln{T^{T, L)) against Tb- Sizes L = 128, 96, 64, 
48, 32, 24, 16 left to right. Same color coding as in Fig. (3] 
Dashed line: extrapolation. Red arrow : exact infinite tem¬ 
perature value. 


with X = Tb we collect together the ThL data points 
for all the sizes L up to x = 0.6 and make standard poly¬ 
nomial fits with 3 or 4 terms y{x) = oq -b aix -I- a 2 X^ 
or y{x) = oo -b aix -b a 2 X^ + a^x^ (in fits with larger 
numbers of terms the fit parameter values become un¬ 
stable). Assuming that each polynomial fit curve ex¬ 
tended to Tf, = 0 is a good approximation to the true 
behavior, each ag is an estimate for the critical limit 
value. The oq values for 3 or 4 parameter fits turn 
out to be similar. In Figs. [HI |T6l [T3 and [18] the data 
and fits are shown for d In x(Tb)/d In Tb, dln(T^(Tb))/dTb 
and 91nx(Tt,)/91n(T^(Tt,)) and dlng(Tb)/dlnTb for both 
Gaussian and bimodal models. The fits are automatic, 
so this procedure is objective and we assume that it is 
optimal for the available data. All the Gaussian extrap- 
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FIG. 13. (Color on line) Gaussian 2D ISG. g{Tb) against 
T(,. Sizes L = 64, 48, 32, 24, 16, 12 top to bottom. Same color 
coding as in Fig. |3] g{T,L) is the Binder cumulant. Line 
slope —3.5. 



FIG. 14. (Color on line) Bimodal 2D ISG. L^g{Tb) against Tb- 
Sizes L = 128, 96, 64, 48, 32, 24, 16, 12 top to bottom. Same 
color coding as in Fig.|3] g{T, L) is the Binder cumulant. Line 
slope —4.2. 
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FIG. 15. (Color on line) The ThL dlnx{Tb)/dlnTb data for 
the Gaussian model from Fig. [SI (lower, black) with the poly¬ 
nomial fit (green), and for the bimodal model from Fig. (6] 
(upper, red) with the polynomial fit (blue). 



dated critical values estimated in this way are close to 
those expected assuming the published exponents, rj = 0 
and v = 3.48(5) [ll|. This implies that the estimated 
bimodal critical values are also close to the true critical 
limits. 


FIG. 16. (Color on line) The ThL i91n(r^(r6)/clr6 data for 
the Gaussian model from Fig. [71 (lower, black) with the poly¬ 
nomial fit (green), and for the bimodal model from Fig. [8] 
(upper, red) with the polynomial fit (blue). 
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